Rods are less fragile than spheres: Structural relaxation in dense liquids composed of 

anisotropic particles 



(N 

o 

(N 
C 

On 



O 

a 

g 
o 

o 



Tianqi Shen 1 , Carl Schreck 1 , Bulbul Chakraborty 2 , Denise E. Freed 3 , and Corey S. O'Hern 4,1 
Department of Physics, Yale University, New Haven, Connecticut 06520-8120, USA 
2 Martin Fisher School of Physics, Brandeis University, 
Mail Stop 057, Waltham, Massachusetts 02454-9110, USA 
3 Schlumberger-Doll Research Center, One Hampshire Street, Cambridge, Massachusetts 02139, USA and 
4 Department of Mechanical Engineering & Materials Science, 
Yale University, New Haven, Connecticut 06520-8286, USA 

We perform extensive molecular dynamics simulations of dense liquids composed of bidisperse 
dimer- and ellipse-shaped particles in 2D that interact via purely-repulsive contact forces. We 
measure the structural relaxation times obtained from the long-time a-decay of the self-part of the 
intermediate scattering function for the translational and rotational degrees of freedom (DOF) as 
a function of packing fraction cj>, temperature T, and aspect ratio a. We are able to collapse the 
packing-fraction and temperature-dependent structural relaxation times for disks, and dimers and 
ellipses over a wide range of a, onto a universal scaling function J r ±{\cj> — (j>o\, T, a), which is similar to 
that employed in previous studies of dense liquids composed of purely repulsive spherical particles in 
3D. T± for both the translational and rotational DOF are characterized by the a-dependent scaling 
exponents /i and S and packing fraction <j>o{a) that signals the crossover in the scaling form J-± from 
hard-particle dynamics to super- Arrhenius behavior for each aspect ratio. We find that the fragility 
at <f>o, m((j>o), decreases monotonically with increasing aspect ratio for both ellipses and dimers. For 
a > a p , where a p is the location of the peak in the packing fraction <f>j at jamming onset, the 
rotational DOF are strongly coupled to the translational DOF and the dynamic scaling exponents 
and 4>o are similar for the rotational and translational DOF. For 1 < a < ct p , the translational 
DOF become frozen at higher temperatures than the rotational DOF, and <f>o for the rotational 
degrees of freedom increases above <j>j. Moreover, the results for the slow dynamics of dense liquids 
composed of dimer- and ellipse-shaped particles are qualitatively the same, despite the fact that zero- 
temperature static packings of dimers are isostatic, while static packings of ellipses are hypostatic. 
Thus, zero-temperature contact counting arguments do not apply to structural relaxation of dense 
liquids of anisotropic particles near the glass transition. 

PACS numbers: 05.10.-l,64.70.Pf,83.80.Fg 



> 

00 



o 

(N 



X 



I. INTRODUCTION 

One of the hallmarks of glassy behavior in dense liquids 
is the rapid increase in the stress and structural relax- 
ation times as the density is increased or the temperature 
is lowered toward the glass transition Several recent 
studies have characterized how the fragility, i. e. the rate 
of increase of the stress and structural relaxation times 
with decreasing temperature near the glass transition [2j , 
depends on the packing fraction in dense liquids of spher- 
ical particles that interact via purely repulsive contact 
potentials in particular, Rcfs. [Ell have shown 

that for dense liquids composed of bidisperse spheres in 
three dimensions (3D) the temperature and packing frac- 
tion dependence of the structural relaxation time r can 
be collapsed onto two master curves, one for <j> > <j>Q and 
another for <f> < 4>q. For 4> < 4>q, the scaling function re- 
duces to hard-sphere behavior r ~ exp [^4/|<^>o — <?^| 5 ] m the 
zero-temperature limit (T — >• 0), where A > is a con- 
stant and S is a scaling exponent that does not depend on 
the form of the repulsive contact potential. For <f> > cf>o, 
the temperature dependence of the structural relaxation 
time is super-Arrhcnius and the fragility increases with 
packing fraction. For 3D bidisperse spheres, (f>o ~ 0.635, 
which is distinct from the jammed packing fraction <f>j for 



both fast and slow packing-generation protocols 

Despite the fact that there have been a number stud- 
ies of the rapidly growing structural and stress relaxation 
times near the glass transition in dense liquids composed 
of anisotropic particles [9l4l3| , there have been few quan- 
titative calculations of the fragility of dense liquids as a 
function of the particle shape. In this article, we inves- 
tigate the slow dynamics of both the translational and 
rotational degrees of freedom in dense amorphous liquids 
composed of bidisperse dimer- and ellipse-shaped parti- 
cles in two dimensions (2D) that interact via purely repul- 
sive contact forces as a function of temperature, packing 
fraction, and aspect ratio. 

We address several important questions: 1) What are 
the forms for the structural relaxation times for the trans- 
lational and rotational degrees of freedom for dimer- 
and ellipse-shaped particles as a function of temperature 
and packing fraction? 2) Can we identify master curves 
for the structural relaxation times for (f> > <po(a) and 
4> < 4>a(a) for each a similar to those found for spher- 
ical particles at a = 1? 3) Are the slow dynamics of 
dense liquids composed of anisotropic, elongated parti- 
cles more or less fragile near 4>o(a) than spherical parti- 
cles? 4) Do differences in microscale geometrical features 
lead to qualitative changes in structural relaxation in sys- 
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terns composed of anisotropic particles? To address these 
questions, we performed extensive molecular dynamics 
simulations of supercooled liquids composed of friction- 
less, purely repulsive dimer- and ellipse-shaped particles 
in 2D over a wide range of temperature T, packing frac- 
tion <j), and aspect ratio a. 




(a) (b) 

FIG. 1: Definition of the aspect ratio a = a/b, where a and 
b are the length of the major and minor axes, for (a) ellipses 
and (b) dimers 



We find four key results. First, we show that the 
structural relaxation times for dimer- and ellipse-shaped 
particles in 2D obey scaling forms similar to that for 
3D systems composed of bidisperse spheres p| Q. In 
particular, we find that there is an aspect-ratio depen- 
dent packing fraction <f>o(a) below which the structural 
relaxation time converges to hard-particle behavior in 
the zero-temperature limit and above which the struc- 
tural relaxation time grows super- Arrheniusly with de- 
creasing temperature. Second, we identify a universal 
scaling form IF±{\4> — 4>o\,T,a) (J r _ for tf> < <fto(a) and 
F+ for (j) > <j>a{a)) that collapses the structural relax- 
ation times for all 2D particle shapes studied, including 
bidisperse disks and both dimer- and ellipse-shaped par- 
ticles over a wide range of a. For small a, the rotational 
degrees of freedom (DOF) are not 'caged' at tempera- 
tures where the translational degrees of freedom become 
frozen, and <j)o(a) for the rotational degrees of freedom 
increases above (j)j(a). For a > a p (the location of the 
peak in 4>j(a) [15[), the rotational degrees of freedom 
are strongly coupled to the translational degrees of free- 
dom, and obey scaling functions with similar scaling ex- 
ponents. Third, the fragility at <po decreases monoton- 
ically with increasing aspect ratio for both dimers and 
ellipses. Finally, we do not find qualitative differences in 
the slow dynamics of dense liquids composed of dimer- 
and ellipse-shaped particles, despite the fact that dimer 
packings at zero temperature are isostatic [l6[ at jam- 
ming onset with on average six contacts per particle in 
2D (z = Zj so = 6) for I < a < 2, whereas ellipse packings 
at zero temperature are hypostatic with z < Zi SO over the 
same range in a [TBj . Naively, one might have expected 
that the structural relaxation times for hypostatic sys- 
tems decay more quickly than those for isostatic systems 
at a given packing fraction \<f> — <pj(a)\ and temperature 
T. However, we show that counting arguments for zero 
temperature packings do not apply for dense liquids com- 
posed of anisotropic particles near the glass transition. 




(a) (b) 

FIG. 2: (a) Schematic of ellipses i and j with orientations 
tii and Uj, respectively, and interparticle separation vector 
fij. 9i is the angle between ti; and the horizontal axis, (b) 
Definition of the contact distance <Jij between ellipses i and 
j, which is obtained by moving ellipse i toward j along 
while fixing their orientations until they are in contact. 



The organization of the remainder of the article is as 
follows. In SecHH we describe the methods we employed 
to measure the structural relaxation times for the trans- 
lational and rotational degrees of freedom for dense liq- 
uids composed of anisotropic particles in 2D. As a test 
of our methods, we show in Sec. [TT] the results for the 
collapse of the structural relaxation times for dense liq- 
uids composed of 2D bidisperse disks onto the scaling 
form F±{\<f>— 0o I, T) used in Refs. [1,0] for 3D bidisperse 
spheres. In Sec. IIIIl we present our novel results for the 
structural relaxation times for the rotational and trans- 
lational degrees of freedom for dimer- and ellipse-shaped 
particles as a function of temperature, packing fraction, 
and aspect ratio, and show that the relaxation time data 
for dimers and ellipses can be collapsed onto a gener- 
alized scaling function F±(\<j> — 4>o\,T,a). In Sec. IIV1 
we summarize our results and identify several important 
open questions for future research studies. The Appendix 
includes additional numerical calculations of the system- 
size dependence of the structural relaxation times and 
nematic order as a function of temperature and packing 
fraction that supplement the results in Sec. IIIIl 

II. METHODS 

We performed molecular dynamics (MD) simulations 
of systems composed of bidisperse, frictionless rigid 
dimer- and ellipse-shaped particles with mass m in 2D 
that interact via purely repulsive contact forces. We 
chose bidisperse systems, where half of the particles are 
large with major axis a; = 1.4 and half are small with ma- 
jor axis a s — 1, to prevent positional and orientational or- 
dering. We varied the aspect ratio (ratio of the major and 
minor axes) a = ai/bi = a s /b s in the range 1 < a < 2.2. 
(See Fig. [TJ) Note that rigid dimers with a > 2 possess 
a small gap between pairs of monomers that comprise a 
given dimer. We employed periodic boundary conditions 
in a square cell with length L and studied system sizes 
in the range N = 64 to 256 particles to assess finite-size 
effects. 
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FIG. 3: Plot of the scaled self-part of the intermediate scat- 
tering function \og[log(Ct/F t (q,t))]/ fit for the translational 
degrees of freedom (symbols) versus \og(t/rt), where Ct is 
the prefactor, fit is the stretching exponent, and 7t is the 
characteristic relaxation time obtained from fits to Eq. [4] for 
bidisperse disks for packing fractions from 0.78 to 0.90 and 
temperature from 3x 10 -6 to 4x 10 -3 . The solid line has slope 
1. (Inset, top) Ft(q,t) versus t (gray lines) with fits (symbols) 
to the stretched exponential form in Eq. [3] in the long-time 
a-decay regime for <f) — 0.83 and temperatures T = 1 X 10~ 3 , 
9 x 10~ 4 , 6 x 10~ 4 , and 5 x 10 4 from left to right. (Inset, 
bottom) The stretching exponent fit from fits of Ft(q,t) to 
Eq. [4] as a function of temperature T at <f> — 0.78 (circles), 
0.83 (triangles), and 0.90 (crosses). 



The particles interact via the pairwise, purely repulsive 
linear spring potential, 

where e is the characteristic energy scale of the repulsive 
interaction, Q(x) is the Heaviside step function, so that 
particles only interact when they overlap, and <Jij is the 
contact distance between particles i and j that in gen- 
eral depends on the particle orientations and separation 
vector fij. For ellipses, fy connects the centers of mass 
of ellipses i and j as shown in Fig. [2j while for dimers, 
fij connects monomers i and j on distinct dimers k and 
I. For dimers, cr^ = (<Tj + aj)/2, where <jj is the di- 
ameter of monomer i. The method for calculating the 
contact distance between two ellipses is discussed in de- 
tail in Ref. [13] • Both the large and small ellipses possess 
the same moment of inertia I e — (ai/a s )mb 2 (l + a 2 )/16, 
while the small dimers have a moment of inertia Ids — 
mb 2 (a — l) 2 /4 and large dimers have a moment of inertia 
Idi = mb 2 (a — l) 2 (ai/2a s ) 2 . We chose e and by/m/e as 
the energy and time units. 

To prepare the systems at a given temperature, we 
first solve Newton's equations of motion using the ve- 




FIG. 4: (Color online) Snapshots of dense liquid configura- 
tions of bidisperse dimer- (left) and ellipse-shaped particles 
(right) with aspect ratio a = 1.8, packing fraction cj> — 0.82, 
and temperature T « 10 -4 . 



locit y V erlet algorithm with a velocity rescaling thermo- 
stat [l8| using a time step At = 0.01. We then equilibrate 
and study the system dynamics without the thermostat 
in the NVE ensemble. To identify the structural relax- 
ation times, we measure the self-part of the intermediate 
scattering function (ISF) for the translational degrees of 
freedom, 

%t) = ^^E eifR(t) ^' (0)1 )' ( 2 ) 

where the wavenumber q = 2n/ai, Ni = N/2 is the num- 
ber of large particles, and (•) indicates an average over 
the orientations of q and time origins. To quantify relax- 
ation of the rotational degrees of freedom, we calculated 

Fr(n ) t) = i^e^W-W°)^ ) (3) 

where 6i is the angle between the orientation il; of particle 
i and the horizontal axis and (•) indicates an average over 
time origins. For most studies, we focused on n = 1, in 
which case relaxation of the ISF occurs when the dimer- 
or ellipse-shaped particles flip by it. 



<t>o € 4>? <P 

FIG. 5: Schematic of the possible packing fractions at jam- 
ming onset 4>J (shaded region) that can be obtained as a 
function of control parameters in packing-generation proto- 
cols {e.g. compression or energy dissipation rates), where 
4> j < <f> a j < 0j tal . The most dilute static packings at (f>j are 
obtained for infinitely fast quenches, while the densest static 
packings at <^>j tal (with either complex unit cells or phase- 
separated packings for bidisperse systems) are obtained for 
infinitely slow quenches. At finite rates, static packings with 
intermediate packing fraction <j>j can be obtained. Note that 
the current studies and those in Refs. [^, 0] indicate that 

0o < 4>j. 
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We measured the structural relaxation times for dense 
liquids composed of bidisperse disks, dimer-, and ellipse- 
shaped particles over a wide range in packing fraction 
from 0.75 to 0.94, temperature from 10~ 2 to 10~ 6 , and 
aspect ratio 1 < a < 2. In this regime, the ISF for the 
translational and rotational degrees of freedom display 
stretched exponential relaxation 

F t (q,t) « Gte-^l^ 1 (4) 

F r (l,t) « C r e- {t/r ^ r (5) 

at long times in the a-decay regime, where Ct and C r 
are prefactors, fit and j3 r are the stretching exponents, 
and Tt and T r are characteristic relaxation times for the 
translational and rotational degrees of freedom, respec- 
tively. In Fig. [3] we show the ISF for the translational 
degrees of freedom for bidisperse disks in 2D with fits 
to the stretched exponential form (Eq. [4j at long times. 
We find that the stretching exponent f3 t decreases with 
increasing packing fraction and decreasing temperature 
(lower inset to Fig. [3]). In Sec. IIII1 we will show results 
from measurements of the decay of the self-part of the in- 
termediate scattering function for both the translational 
and rotational degrees of freedom for dimer- and ellipse- 
shaped particles as a function of (f>, T, and a. Prior to all 
measurements, we run MD simulations of the system for 
10 times the longest relaxation time r to reach metastablc 
equilibrium. 
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FIG. 6: (Color online) (a) Structural relaxation time r t \fT 
versus 1/T for dense liquids of bidisperse disks in 2D for 23 
different packing fractions, <j> = 0.94, 0.92, 0.90, 0.88, 0.86, 
0.85, 0.845, 0.84, 0.835, 0.83, 0.825, 0.82, 0.815, 0.81, 0.805, 
0.80, 0.795, 0.79, 0.785, 0.78, 0.77, 0.76, and 0.75 from left 
to right, (b) Collapse of the data in (a) using Eq. [7] with 
4>o — 0.831, fi — 1.25, and 8 = 1.95. The dashed line has 
slope fiS/2. Inset: Scaling relation between the total potential 
energy and temperature V = VqT^ at <j> = <f>o ~ 0.83, where 
Vo — 2.0 and fi = 1.25. The dashed line has slope 1.25. 



TABLE I: Dynamic scaling exponents fi and 8 obtained from 
fitting the structural relaxation time Tt (<j), T) for the trans- 
lational degrees of freedom to the form in Eq. [7] fragility 
m{4>o) = [i8j2 at cf>o, and the packing fractions (f>o, (f>j, and 
4>j for dense liquids composed of bidisperse, purely repulsive 
disks in 2D and spheres in 3D [5J. 

We also monitor the positional and orientational or- 
dering of the systems as a function of decreasing temper- 
ature by measuring the bond orientational order [l9j and 
nematic order parameter 

P 2 = (2cos 2 (9 l -e)-l), (6) 

where 8 is the average particle orientation in a given con- 
figuration and the angle brackets indicate an average over 
independent configurations. As shown in the snapshots 
of the configurations in Fig. 31 we do not find signifi- 
cant positional or orientational ordering in dense liquids 
of bidisperse dimer- and ellipse-shaped particles. 

In previous studies of slow dynamics in dense liquids, 
Witten and Berthier Q identified a dynamical scaling 
relation, 

rW ,,T)VT~exp^— W F ± ^ jj, (7) 



which collapsed the structural relaxation time t(</>, T) 
from the self-part of the intermediate scattering function 
for bidisperse, purely repulsive spheres in 3D over a wide 
range of temperature and packing fraction. The scaling 
function F± assumes two different forms above and below 
4>q ~ 0.635, and is only a function of \(j>o~4>\ 2 ^ /T , where 
\x is an exponent that controls the scaling of the poten- 
tial energy with temperature, V ~ VqT^, near cj>g. These 
studies showed that r t displays hard-sphere dynamics in 
the zero-temperature limit for <j) < <j>$, but r t displays 
super-Arrhenius behavior with decreasing temperature 
for (j) > 4> . For 4> -t cf> , x = \<f> - <f>\ 2 fi i /T < 1 and 
F±(x) ~ x 5 ^' 2 , which implies that r - exp^/T^/ 2 ) 
and the fragility m(0o) = ^5/2 at (j)o for T — > is con- 
trolled by the scaling exponents /i and S. The scaling 
exponent S does not depend of the form of the purely 
repulsive contact potential. The results for the scaling 
exponents \x and S, (f>o, and fragility m((f>o) from Ref. [5j 
are shown in Table U for dense liquids composed of bidis- 
perse, purely repulsive spheres in 3D. 

In Ref. Berthier and Witten also compared the 
critical packing fraction <f>Q ~ 0.635 (which signals the 



FIG. 7: Ensemble averaged packing fraction <j>j at jamming 
onset versus aspect ratio a for N = 480 bidisperse, frictionless 
dimer- (squares) and ellipse-shaped (circles) particles. 



crossover from hard-sphere dynamics to super- Arrhenius 
temperature dependence of r) to the packing fraction (f>j 
at jamming onset. Numerous studies have shown that 
4>j for frictionless, purely repulsive spherical particles 
depends on control parameters in the jammed packing- 
generation protocol (ej7. the compression and energy 
dissipation rates) [8, 20-22] as shown schematically in 
Fig. [5] In the limit of infinitely fast rates, the most dilute 
and disordered mechanically stable packings at 4>j are ob- 
tained. In contrast, in the limit of infinitely slow rates, 
the densest crystalline packings are obtained. At finite 
rates, intermediate values of the packing fraction at jam- 
ming onset and positional and other types of order can 
be obtained. Previous studies have shown for bidisperse 
systems that disordered mechanically stable packings at 
jamming onset can exist over a finite range of packing 
fractions (i.e. from cf>j to cf> s j in Fig. [SJ with only small 
and subtle changes in positional order. For 3D bidisperse 
spheres, 0.648 f| and <j> s j w 0.662 @. The scaling 

analyses in Refs. [5j, |6| for the structural relaxation times 
for 3D bidisperse systems suggest that fa < <f)j. 

Recent experimental [H, [24| and computational [H, 
[25j studies have shown that the average jammed pack- 
ing fraction (in the limit of fast quenching rates) (j)j{a) 
for static packings of anisotropic frictionless particles in- 
creases with aspect ratio a from • random close pack- 
ing' [26J at a = 1, reaches a peak near a p ~ 1.5 (where 
a p depends on the specific particle shape and spatial di- 
mension), and decreases continuously beyond a p . (See 
Fig. [7] for a plot of the jammed packing fraction <fij(a) 
for ellipses and dimers as a function of aspect ratio in 
two dimensions (2D).) In these studies, the jammed pack- 
ings were created using 'fast' packing-generation proto- 
cols, and thus these packings do not possess significant 
positional or orientational order [Tt], H3, [28| . 
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FIG. 8: (Color online) Self-part of the intermediate scattering 
function (ISF) for the (a) translational and (b) rotational (n = 
1) degrees of freedom for dimers at aspect ratio a = 1.6 and 
packing fraction (f> = 0.86 over a range of temperatures T = 
5 x 10" 3 , 3.3 x 10" 3 , 2.2 x 10" 3 , 1.4 x 10" 3 , 9.3 x 10" 4 , 
6.2 x 10~ 4 , 4 x 10~ 4 , 2.7 x 10~ 4 , 1.8 x 10~ 4 and 1.2 x 10~ 4 
from left to right. The dashed horizontal line in (a) indicates 
the approximate value of the plateau in the ISF. The symbols 
in (a) and (b) show fits of the ISF to stretched exponential 
behavior (Eqs. |4]and[5)l for long times. 



III. SCALING ANALYSIS FOR THE 
STRUCTURAL RELAXATION TIMES 

In this section, we describe our measurements of the 
structural relaxation time r from the self-part of the in- 
termediate scattering function for the translational and 
rotational degrees of freedom for bidisperse disks, dimer- 
and ellipse-shaped particles. We will investigate the ex- 
tent to which the scaling relation in Eq. [JJ holds for 
2D systems as a function of aspect ratio and determine 
the variation with aspect ratio of the fragility at fa (a) 
that signals the crossover from hard-particle dynamics 
to super-Arrhenius temperature dependence. We find no 
significant positional or nematic ordering with decreasing 
temperature over the full range of packing fraction and 
aspect ratio. (See Fig. [16] that shows the nematic order 
parameter P2 < 1/y/N for a = 1.8 ellipses over a wide 
range of packing fraction and temperature.) 

We first focus on the slow dynamics of dense liquids 
composed of bidisperse disks in 2D near jamming to com- 
pare our results with those obtained for dense 3D liquids 
of bidisperse spheres [j| 0. The structural relaxation 
time Tt obtained from the self-part of the intermediate 
scattering function is plotted versus the inverse temper- 
ature 1/T in Fig. [S] (a) over a wide range in packing 
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FIG. 9: (Color online) Self-part of the intermediate scattering 
function for the (a) translational and (b) rotational degrees 
of freedom (n = 1) for ellipses at aspect ratio a = 1.6 and 
packing fraction (j> — 0.87 over a range of temperature T — 
4x 10~ 3 , 3x 10~ 3 , 1.8xl0~ 3 , 1.4xl0 -3 , 1 x 1CT 3 , 6.2 x 10 -4 , 
3.7 x 1(T 4 , 2.3 x 10~ 4 , 1.7 x 10~ 4 , 1.2 x 1(T 4 , and 1.0 x 10~ 4 . 
The dashed horizontal line in (a) indicates the approximate 
value of the plateau in the ISF. The symbols in (a) and (b) 
show fits of the ISF to stretched exponential behavior (Eqs. [4] 
and[5J) for long times. 
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FIG. 10: Plot of the scaled self-part of the intermediate scat- 
tering function log [log (C/F(q, t))]//3 for the translational and 
rotational degrees of freedom (symbols) versus \og(t/r), where 
C is the prefactor, (3 is the stretching exponent, and r is 
the characteristic relaxation time obtained from fits to Eqs. [4] 
and[5]for ellipses and dimers at packing fraction <f> = 0.87 over 
a range of temperature from 2 x 10 -3 to 1 x 10 -4 . The solid 
line has slope 1. (Inset, top) F t (q,t) versus t (gray lines) with 
fits (symbols) to the stretched exponential form in Eq. [4] in 
the long-time a-decay regime for ellipses at <f> = 0.87 and tem- 
peratures T = 1.8 X 10" 3 , 1.4 x 10" 3 , 1.0 x 10" 3 , 6.2 x 10" 4 , 
3.7 x 10" 4 , 2.3 x 10" 4 , and 1.0 x 10~ 4 from left to right. (In- 
set, bottom) The stretching exponent fit from fits of Ft(q,t) 
to Eq. [4] as a function of temperature T at <f> = 0.87 (circles) 
for ellipses. 



fraction from <f> = 0.75 to 0.94. One can clearly identify 
a change in the form of the structural relaxation time as 
the packing fraction increases above <fio ~ 0.82 — 0.84. 
In Fig. [6] (b), we collapse the data in (a) using the scal- 
ing form in Eq. [7] by plotting \<po — 4>\ s \og(r t VT) ver- 
sus \4>o ~ <j>\ 2 '^/T. The exponent fi is obtained from the 
power-law scaling of the total potential energy V with 
temperature for <j> « <fio @, 0| . In the inset to Fig. [5] (b) , 
we show that V ~ T M over several orders of magnitude 
in temperature with an exponent fi = 1.25, which is sim- 
ilar to the corresponding value for 3D bidisperse spheres, 
/i3D = 1.3 jj}. The scaling form for the structural relax- 
ation time possesses two distinct branches: F + for <f) > <fio 
and F_ for <j> < 4> . To obtain Fig. [6] (b), we chose the 
combination of parameters <j) = 0.831, fJ,(<po) — 1-25, 
and 8 = 1.19 that yielded the tightest collapse of the Tt 
data in Fig. [6] (a). For 2D bidisperse systems, we find 
that the fragility at </> , m(0 o ) = 1.2, which is slightly 
lower than the comparable value, m(4>o) — 1.4, for 3D 
bidisperse spheres. (See Table HI) 

In Figs. H] and HO we show the ISF for the translational 
and rotational degrees of freedom for dimers and ellipses, 
respectively, over a range of temperatures at comparable 
values of the packing fraction <fi — 4>j(a) for aspect ra- 
tio a = 1.6 above a p at which the peak in <j)j(a) occurs 



(Fig. [7]). For both dimers and ellipses, the decay of the 
ISF for the translational and rotational degrees of free- 
dom at long-times can be described using a stretched 
exponential form as shown in Fig 1101 For a > a p the 
rotational degrees of freedom are strongly coupled to the 
translational degrees of freedom and r r ~ r t . For both 
dimers and ellipses, the 'plateau' value of the ISF is larger 
for the rotational degrees of freedom compared to that for 
the translational degrees of freedom. The only significant 
difference between the intermediate scattering functions 
for dimers and ellipses in Figs.[S]and|n]is that the plateau 
value for the translational degrees of freedom is smaller 
(F t ~ 0.6 - 0.7) compared to that for dimers (F t ~ 0.95), 
which indicates that the 'cage' size for ellipses is larger 
than that for dimers. 

We use a similar procedure to that employed for bidis- 
perse disks to collapse the data for the structural relax- 
ation times t> and r t for dimer- and ellipse-shaped par- 
ticles using the scaling form in Eq. [7J Examples of the 
collapse of the structural relaxation times for the trans- 
lational degrees of freedom for dimers and ellipses are 
shown in Figs. [TT] and [T^l respectively, at aspect ratio 
a = 1.8. Similar collapse for r r is found for both dimer- 
and ellipse-shaped particles. Also, the quality of the data 
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FIG. 11: (Color online) (a) Structural relaxation time for the 
translational degrees of freedom r t VT as a function of inverse 
temperature 1/T on a log-log scale for dimers at aspect ratio 
a = 1.8 over a range of packing fractions <f> — 0.95, 0.93, 0.91, 
0.90, 0.89, 0.88, 0.875, 0.87, 0.865, 0.86, 0.855, 0.85, 0.845, 
0.84, 0.835, 0.83, 0.825, 0.82, 0.815, 0.81, 0.80, 0.79, and 0.78 
from left to right, (b) Collapse of the data in (a) using Eq. [7] 
with 0o = 0.852, n = 1.29, and 5 = 0.85. The dashed line has 
slope n8/2. 



collapse is the same for dimers and ellipses. We plot 
the characteristic packing fraction <fio (a) that signals the 
change in scaling form for t(</>, T), the scaling exponents 
/i and <5, and the fragility m(<fio) for the rotational and 
translational degrees of freedom as a function of aspect 
ratio a for dimers and ellipses in Figs.[T3land ll4[ respec- 
tively. 

For dimers, the fragility m(0o) for the translational 
degrees of freedom decreases by a factor of ~ 2.5 with in- 
creasing aspect ratio from m((f>o) « 1.2 at a = 1 to « 0.5 
at a = 1.8. The decrease in fragility is caused by the de- 
crease in the scaling exponent 5 (and thus insensitive to 
the form of the purely repulsive contact forces) because 
/.t 1.25 is nearly constant over this range in aspect ra- 
tio. The fragility for the rotational degrees of freedom 
decreases by a similar factor to that for the translational 
degrees of freedom for a > 1.2. In addition, in Fig.H31(a) 
we show that the a-dependence of </>o for the translational 
degrees of freedom matches that for </>j, i.e. nonmono- 
tonic dependence on a with a peak near a p « 1.4, and 
(f>o satisfies 4>o(&) < 4>j(pt)- For a > a p , 4>o is the same 
for the translational and rotational degrees of freedom. 
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FIG. 12: (Color online) (a) Structural relaxation time for the 
translational degrees of freedom r t VT as a function of inverse 
temperature 1/T on a log-log scale for ellipses at aspect ratio 
a — 1.8 over a range of packing fractions <f> = 0.98, 0.96, 0.94, 
0.92, 0.91, 0.90, 0.89, 0.88, 0.87, 0.86, 0.855, 0.85, 0.84, 0.83, 
0.82, and 0.81 from left to right, (b) Collapse of the data in 
(a) using Eq. [TJwith (j> = 0.884, ^ = 1.26, and 8 = 0.90. The 
dashed line has slope fj.S/2. 



However, for a < a p , (f>Q for the rotational degrees of 
freedom increases with decreasing a and becomes larger 
than 4>q for the translational degrees of freedom and <f>j. 
This strong increase in <po for the rotational degrees of 
freedom at small a indicates the existence of a 'rotator', 
amorphous, solid-like regime [29[ at low temperatures, 
where the translational degrees of freedom are frozen, 
but the rotational degrees of freedom are liquid-like. We 
show that these results are nearly independent of system 
size by comparing the scaling collapse of the structural 
relaxation times for the translational degrees of freedom 
for systems composed of N = 64 and 200 dimers at as- 
pect ratio a — 1.3 in Fig. [171 

We find quantitatively similar behavior for the depen- 
dence of the scaling exponents \x and d, fragility m(0o), 
and packing fraction 4>q on aspect ratio for ellipses and 
dimers as shown in Fig. 1141 A minor difference between 
dimers and ellipses is the weaker decrease in </>o and <pj 
with increasing aspect ratio a. Again, c/jq for the rota- 
tional degrees of freedom increases above 4>o for the trans- 
lational degrees of freedom and <fij for a < a pi which 
indicates the onset of the 'rotator' regime for low aspect 
ratios. In Fig. [T5] we show that we are able to collapse 
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FIG. 13: (Color online) (a) The characteristic packing frac- 
tion (j>o that signals the crossover from hard-particle dynam- 
ics to super- Arrhenius temperature dependence for the struc- 
tural relaxation times for the translational (circles) and rota- 
tional (downward triangles) degrees of freedom and the pack- 
ing fraction <j>J at jamming onset (dashed line) for dimers as 
a function of aspect ratio a. (b) The scaling exponents fi 
(circles, solid line) and 8 for translational (squares, solid line) 
and rotational (diamonds, dashed line) degrees of freedom, 
and fragility m((j>o) for the translational (asterisks, solid line) 
and rotational (triangles, dashed line) degrees of freedom for 
dimers versus aspect ratio a. 
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FIG. 14: (Color online) (a) The characteristic packing frac- 
tion c/>o that signals the crossover from hard-particle dynamics 
to super- Arrhenius temperature dependence for the structural 
relaxation times for the translational (circles) and rotational 
(downward triangles) degrees of freedom and the packing frac- 
tion cf)j at jamming onset (dashed line) for ellipses as a func- 
tion of aspect ratio a. (b) The scaling exponents /i (circles, 
solid line) and <5 for translational (squares, solid line) and 
rotational (diamonds, dashed line) degrees of freedom, and 
fragility m(4>o) for the translational (asterisks, solid line) and 
rotational (triangles, dashed line) degrees of freedom for el- 
lipses versus aspect ratio a. 



all of the structural relaxation time data for the trans- 
lational and rotational degrees of freedom for bidisperse 
dimers and ellipses using the generalized scaling form 

F±{x) = c(a)[F ± (x)} 1/m ^ o(a) \ (8) 

where c(a) is an a-dependent prefactor that is differ- 
ent for the translational and rotational degrees of free- 
dom and x = \<j> - (/> | 2/Al /T. Thus, we find that the 
scaling behavior of the structural relaxation times for 
dimers and ellipses is nearly identical over a wide range 
of aspect ratios, despite the fact that static packings of 
dimer-shaped particles are isostatic while static packings 
of ellipse-shaped particles are hypostatic at zero temper- 
ature. 



IV. CONCLUSIONS 

We performed extensive molecular dynamics simula- 
tions of the slow dynamics in dense liquids composed 



of bidisperse, purely repulsive dimer- and ellipse-shaped 
particles in 2D. We showed that a similar scaling form 
to that identified for structural relaxation for spherical 
particles @ is able to collapse the packing fraction 
and temperature-dependent structural relaxation times 
for the rotational and translational degrees of freedom 
for dimer- and ellipse-shaped particles over a wide range 
of aspect ratios a (Fig. [TS]). Thus, the dynamical criti- 
cal point at T — and (f) = <po studied in Refs. 0, @ for 
dense liquids composed of purely repulsive spheres can be 
generalized to T = 0, <j> = </>o ; an d a — 1, and this dynam- 
ical critical point controls the temperature, packing frac- 
tion, and aspect-ratio-dependent fragility of dense liquids 
composed of dimer- and ellipse-shaped particles. Fur- 
thermore, we find qualitatively and quantitatively simi- 
lar results for the dynamical scaling exponents, packing 
fraction (f>o that signals the crossover from hard-particle 
dynamics to super-Arrhenius temperature dependence, 
and fragility m((f>o) for dimer- and ellipse-shaped parti- 
cles. In particular, for dimers and ellipses the fragility 
at 4>q decreases monotonically with increasing aspect ra- 
tio. Thus, the microscale differences in shape between 
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FIG. 15: The scaling function J-±(x) defined in Eq. [8] versus 
x = \<j> — 0o| 2/ ' M /T ', that collapses the structural relaxation 
times for the translational (Figs. [TT] (b) and 1 121 (b)) and rota- 
tional degrees of freedom for ellipses and dimers over a wide 
range of aspect ratios. The dashed line has slope 1. Insets: 
c(a) for the translational (solid) and rotational (dashed) de- 
grees of freedom for ellipses (top) and dimers (bottom). 



dimers and ellipses do not give rise to important differ- 
ences in the slow dynamics of structural relaxation in 
dense supercooled liquids of anisotropic particles. In con- 
trast, static packings of dimers at zero temperature are 
isostatic, while static packings of ellipses are hypostatic, 
which leads to significant differences in the density of 
vibrational modes and static shear modulus at zero tem- 
perature fl7l.l27l|. 

These findings suggest several important future stud- 
ies. First, we will investigate possible causes of the de- 
crease in fragility of dense liquids composed of elongated 
particles with increasing aspect ratio, including an in- 
crease in the number of inherent structures [30j and vi- 
brational entropy [l[ with increasing a. Second, is there 
a nonequilibrium, finite-temperature regime where mi- 
croscale differences between anisotropic particle shapes 
can cause significant differences in the time-dependent 
structural and mechanical response of the system? To 
address this question, we will subject jammed pack- 
ings of dimers and ellipses to weak thermal fluctuations 
(much below the temperatures studied in the present 
manuscript) and measure the time-dependent structural 
and stress relaxation. 
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Appendix A: Supplementary Results 

In this section, we provide additional numerical calcu- 
lations that supplement the results presented in Sec. IIIII 
In Fig. [TH1 we show the average nematic order Pi (Eq. [5]) 
for dense liquids composed of ellipses with a = 1.8 ver- 
sus l/T over a wide range of packing fractions above and 
below tfio. Pi < 1/vN for all systems studied, which 
indicates that the ellipse orientations are disordered. We 
obtain similar results to those in Fig. [T5] for dense liq- 
uids of bidisperse dimers over the range of T, <f>, and a 
considered. 
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FIG. 16: Nematic order parameter P2 versus log(l/T") for 
dense liquids composed of N — 200 ellipses at aspect ra- 
tio a = 1.8 averaged over 100 configurations for a range of 
packing fractions <j> = 0.81 (open circles), 0.82 (squares), 0.83 
(leftward triangles), 0.84 (filled circles), 0.85 (pluses), 0.86 
(asterisks), 0.87 (diamonds), 0.88 (downward triangles), and 
0.90 (upward triangles). The dashed horizontal line indicates 
P 2 = 1/VN. 



In Fig. [T71 we investigate the effect of system size on 
the structural relaxation times. We show that the scaled 
structural relaxation times for the translational degrees 
of freedom for dense liquids composed of dimers for two 
system sizes (N = 64 and 200) are nearly identical. Sim- 
ilar results to those in Fig. [T7] for dimers hold over the 
full range of aspect ratio considered. We also find that 
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there are negligible system-size effects for the structural 
relaxation times for ellipses with N > 64 for 1 < a < 2. 

10 1 r , , , -, , , 




FIG. 17: Scaled structural relaxation times for the transla- 
tional degrees of freedom \</>o — <j>\ logrtVT versus the scaled 
temperature \(f>o — 0| 2 ' M /T on a log- log scale for dimers at 
a — 1.3 for N = 64 (squares) and 200 (circles) particles. The 
system-size dependence of the scaling parameters is weak: 
0o = 0.878 (0.876); fi = 1.23 (1.25); and S = 1.1 (1.1) for 
N — 256 (64) particles. The dashed line has slope jj,8/2 for 
TV = 256. 
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